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PREFACE 


This  study  was  conducted  under  DA  Project  4A161 102B52C,  Task  B,  Work 
Unit  012,  "Electronic  Image  Analysis  for  Feature  Extraction." 
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Center  for  Physical  Sciences;  and  Mr.  M.  Crowell,  Jr.,  Director,  Research 
Institute. 

COL  Edward  K.  Wintz,  CE,  was  Commander  and  Director  and  Mr.  Robert  P. 
Macchia  was  Technical  Director  of  the  Engineer  Topographic  Laboratories 
during  the  study  period. 
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RADAR  BRIDGE  PATTERNS  EXTRACTION  AND  RECOGNITION 


INTRODUCTION 

The  objective  of  this  study  is  to  devise  a  feature  extractor  for 
extracting  various  bridge  patterns  from  radar  imagery.  This  extractor  can 
be  used  in  conjunction  with  other  feature  extractors  developed  earlier, 
for  example,  the  feature  extractor  using  Walsh  transforms  for  recognizing 
road  intersections,  line  roads,  and  rectangular  objects1  and  the  feature 
extractor  employing  image  statistics  for  classifying  a  set  of  imagery  such 
as  city,  water,  field  and  forest. ^ 

The  technique  presented  makes  use  of  a  combination  of  the  Hough  transform, 
the  Sobel  edge  operator,  image  thresholding  methods,  and  decision  logics 
to  detect  and  recognize  bridge  candidate  patterns  that  are,  in  general, 
characterized  by  a  pair  of  long  lines  that  are  relatively  close  to  each 
other. 

The  hardware  system  consists  of  a  32-element  by  32-element  solid-state 
array  to  convert  the  bridge  patterns  into  electronic  signals,  a 
minicomputer  to  process  the  electronic  signals  from  the  array  using  a 
controlled  software,  and  a  computer-controlled,  two-dimensional  stage  as 
the  image  holder.  The  majority  of  the  bridge  detection  and  recognition 
routines  have  been  implemented  as  software  for  the  hardware  system 
described. 

The  description  of  the  system  operation  is  followed  by  an  explanation  of 
the  strategies  of  bridge  detection  and  recognition  and  bridge  orientation 
angle  determination.  The  system  test  results  are  discussed  and 

conclusions  presented. 


SYSTEM  DESCRIPTION 

The  block  diagram  of  the  system  is  shown  in  figure  1.  A  9-  by  9-inch 
glass  plate  mounted  with  various  2-  by  2-inch  cut  radar  imagery  is 
illuminated  by  a  light  source,  and  a  section  of  the  image  containing  a 


^P.F.  Chen,  F.W.  Rohde,  and  W.W.  Seemuller,  Classification  of  Cartographic 
Features  Through  Walsh  Transforms,  U.S.  Army  Engineer  Topographic 
Laboratories,  Fort  Belvoir,  Va.,  ETL-0290,  April  1982,  AD-116  731. 

2 

P.F.  Chen,  Preliminary  Radar  Feature  Extraction  and  Recognition  Using 
Texture  Measurement ,  U.S.  Army  Engineer  Topographic  Laboratories,  Fort 
Belvoir,  Va.,  Research  note  in  progress. 


bridge  pattern  is  projected  onto  the  surface  of  a  Reticon  32-element  by 
32-element  solid-state  array  through  an  imaging  lens.  The  array  converts 
the  optical  energy  of  the  bridge  pattern  into  a  video  signal.  The  video 
signal  is  quantized  into  10  bits  of  digital  signals  and  sent  to  the 
Hewlett-Packard  2117F  minicomputer  for  processing.  The  computer  first 
takes  in  a  frame  of  32  by  32  pixels  of  the  digitized  signal.  Each  pixel 
is  sequentially  compared  to  an  automatically  set  threshold  value.  It  is 
set  to  a  "1"  if  it  is  greater  than  the  threshold  value  and  to  a  "0" 
otherwise. 


FIGURE  1.  System  Block  Diagram. 

The  edges  of  the  signal  after  thresholding  are  enhanced  using  the  Sobel 
edge  operator.  The  Hough  transform  of  the  edge-enhanced  signal  is  then 
computed.  The  transformed  matrix  has  a  dimension  of  64  by  32  and  displays 
unique  characteristics  of  a  bridge.  A  recognition  routine  follows  to 
screen  out  this  unique  characteristic  and  classifies  the  image  as  a  bridge 
with  its  orientation  angle  against  the  X-axis  of  the  measurement  window 
(or  the  solid-state  array).  In  case  the  unique  characteristic  is  missing, 
the  image  is  identified  as  "not  a  bridge."  The  recognition  result  is 
indicated  on  the  CRT  console.  At  the  end  of  classification,  the  computer 
sends  a  signal  to  the  translational  stage  controllers  to  move  the  stages 
in  the  predetermined  x  and  y  positions,  and  a  new  section  of  image  is 
projected  onto  the  surface  of  the  solid-state  array.  The  procedure 
described  repeats  until  all  preselected  images  are  classified. 


BRIDGE  DETECTION  AND  RECOGNITION 


As  mentioned  earlier,  a  candidate  bridge  pattern,  in  general,  consists  of 
a  pair  of  long  lines  that  are  relatively  close  together.  Thus,  each 

possible  bridge  pattern  is  subjected  to  a  set  of  tests  that  determine  the 
pattern  to  be  classified  as  either  "a  bridge"  or  "not  a  bridge." 

The  process  involves  at  least  the  following  four  steps: 

1.  Preprocessing  and  thresholding  to  eliminate  or  reduce  background 
and  unwanted  noise. 

2.  Edge  detection  to  produce  a  pair  of  long,  parallel  lines. 

3.  Transformation  of  the  long,  parallel  lines  into  a  representation 
more  easily  identified  by  the  computer. 

4.  Test  to  determine  pattern  as  "a  bridge"  or  "not  a  bridge." 

To  eliminate  or  reduce  the  background  and  unwanted  noise,  the  incoming 

pixels  were  compared  to  an  automatically  set  threshold  value.  The  pixels 
that  were  greater  than  the  threshold  value  were  set  equal  to  "1"  and  those 
that  were  less  than  or  equal  to  the  threshold  value  were  set  equal  to 

"0."  The  threshold  value  was  set  equal  to  a  half  of  the  maximum  pixel 
value  within  a  given  frame  of  pixels.  With  this  method  the  unwanted 

background  noise  was  eliminated  in  most  of  the  cases. 

The  edge  enhancement  was  accomplished  by  using  a  Sobel  operator  to  produce 
a  pair  of  long,  parallel  edge  lines  of  the  bridge  pattern.  The  Sobel 
operator,  according  to  our  past  experience,  has  proved  to  be  the  best  for 
generating  clear  and  faithful  image  edges. 

Detection  of  the  pair  of  long,  parallel  lines  produced  by  the  Sobel 
operation  in  the  spatial  domain  was  relatively  complicated.  Therefore, 
the  lines  were  transformed  into  the  p  and  0  domain  using  the  Hough 
transform.  The  Hough  transform  maps  a  line  in  Cartesian  coordinate  space 
(x  and  y  spatial  domain)  into  a  point  in  polar  coordinate  space 
(Pand  0).  A  family  of  lines  with  different  slopes,  passing  througha 
common  point  in  the  x-y  domain,  transforms  into  a  connected  set  of  P-0 
points.  For  our  case,  each  edge  point  produced  by  the  Sobel  operator  in 


■> 

W.K.  Pratt,  Digital  Image  Processing,  New  York:  John  Wiley  and  sons, 
Inc.,  1978. 

^W.K.  Pratt,  Digital  Image  Processing,  New  York:  John  Wiley  &  Sons,  Inc., 
1978. 
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the  x-y  domain  was  transformed  to  a  curve  in  the  P-0  domain,  which  was 

quantized  into  cells.  If  an  element  of  a  curve  fell  in  a  cell,  that 

particular  cell  was  increased  by  one  count.  After  all  edge  points  were 

transformed,  the  P-0  cells  were  examined.  Large  cell  counts  correspond 
to  colinear  edge  points  that  may  be  fitted  by  a  straight  line  with  the 
appropriate  (P,0)  parameters.  The  size  of  the  transform  matrix  was 
selected  to  be  64  by  32  to  yield  enough  resolution  for  determining  the 

orientation  of  bridge  patterns.  The  maximum  available  matrix  size  of  the 
Hough  transform  is  directly  related  to  the  size  of  the  measurement  window 
(input  sampling  array). 

To  determine  the  orientation  angle  of  the  bridge  patterns  against  the  X- 
axis  of  the  measurement  window,  the  10  largest  numbers  were  sorted  in  a 
decreasing  magnitude  order  together  with  their  indices  (the  T's  and  J's 
specifying  the  location  of  the  transform  element  holding  the  sorted 
numbers)  from  the  transform  matrix.  These  10  numbers  were  compared  to  a 
second  threshold  value,  "ITT,"  which  was  automatically  set  equal  to  0.75 
of  the  maximum  element  value  of  the  Hough  transform.  Only  the  numbers 
greater  than  "itt"  were  screened  and  saved  from  recognition  and  detection 
of  bridge  patterns  and  their  orientation  angles. 

The  determination  of  the  bridge  orientation  is  as  follows: 

(a)  First,  the  row  number,  J,  of  the  largest  screened  number  is 
found.  If  there  is  at  least  one  other  screened  number  in  the  same  row,  or 
in  other  words,  having  the  same  J,  then  that  J  is  used  for  the  computation 
of  the  bridge  orientation  angle. 

(b)  If  the  condition  (a)  is  not  met,  then  look  for  a  new  row  with  two 
or  more  screened  numbers.  Use  that  row  number,  J,  for  the  computation  of 
the  bridge  orientation  angle. 

(c)  If  neither  condition  (a)  nor  condition  (b)  is  present  then  just 
use  the  row  number,  J,  in  which  the  maximum  screened  number  resides  for 
the  computation  of  the  bridge  orientation  angle. 


-  sjp 


(d)  The  formulas  for  the  computation  of  bridge  angle,  Rj ,  are 
Rj  =>  22.5  ( 32-J)  degrees;  if  20  £  J  £  52  (1) 

R1  *  180  -  2.5  J  degrees;  if  1  _<  J  <  20  (2) 

The  test  for  deciding  whether  a  pattern  belongs  to  the  category  of  a 
bridge  is  as  follows:  Since  a  bridge  pattern  always  displays  some  width 
in  its  appearance  in  radar  imagery,  each  pixel,  after  first  thresholding 
in  the  measurement  window,  is  tested  for  coexistence  with  its  neighbors. 
As  indicated  before,  the  measurement  window  after  first  thresholding 
contains  pixels  of  only  two  values,  either  a  "1"  or  a  "0".  Thus,  if  a 
pixel  itself  is  a  "1"  and  its  left-hand  neighbor  is  also  a  "1",  a  row 
count  will  be  made.  Similarly,  if  a  pixel  is  a  "1"  and  the  pixel  just 
below  it  is  also  a  "1",  then  a  column  count  will  be  made.  The  row  and 
column  counts  are  examined  after  the  tests  are  made  for  all  pixels  in  the 
measurement  window.  A  flag  is  set  if  any  one  of  these  two  counts  is 
greater  than  eight.  This  flag  is  designated  as  "ANS."  This  flag, 
together  with  the  existence  of  one  of  the  first  three  conditions  stated  in 
the  previous  paragraph,  constitutes  the  absolute  requirement  for 
recognizing  and  detecting  bridge  patterns  and  their  orientation  angles. 


SYSTEM  TEST  RESULTS 

A  set  of  high  quality,  scale  of  1  to  100,000,  X-band  synthetic  aperture 
radar  imagery  from  the  Hudson  River,  New  York  State,  area  consisting  of 
bridge  images  (transportation  DIMS  category  number  503,  FIC  260)  was  used 
for  this  experimentation. 

Figures  2(a)  to  9(a)  show  the  line  printer  output  of  bridge  patterns  with 
different  orientation  angles.  Each  pattern  is  printed  in  16  gray  shades 
by  line  printing.  Figures  2(b)  through  9(b)  show  the  recognition  results 
of  bridges  at  various  orientation  angles.  The  Rough  tranpform  for  these 
bridge  patterns  is  illustrated  in  figures  2(c)  to  9(c).  It  is  seen  that 
each  bridge  pattern  was  recognized  properly  and  its  orientation  angle 
identified  correctly.  Figures  10(a),  10(b),  and  10(c)  demonstrate  the 
input  image,  the  recognition  result,  and  the  Hough  transform  when  the 
measurement  window  sampled  an  entirely  black  area.  As  expected,  the 
system  recognized  this  image  as  "not  a  bridge."  The  resolution  of  the 
bridge  detection  with  present  system  setup  is  approximately  22.5 
degrees.  This  resolution  can  be  increased  by  using  a  larger  size 
measurement  window.  Figures  11  and  12  illustrate  the  bridge-recognition 
results  for  the  case  of  using  64-  by  64-  pixels  measurement  window.  As 
expected,  the  bridge-detection  resolution  was  increased  to  approximately 
10  degrees. 
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It  was  discovered  that  the  bridge  orientation  angle  determination  was  less 
accurate  for  bridge  patterns  oriented  at  45  degrees  and  135  degrees  than 
for  other  orientations.  This  is  due  to  the  quantization  error  of  the 
small-size  measurement  window  (32  by  32  pixels  for  our  case)  that  distorts 
the  straight  bridge-edge  lines  into  looking  like  staircase  patterns  at  45 
degrees  and  135  degrees,  especially  when  the  patterns  shift  their  position 
from  one  corner  of  the  measurement  window  to  the  other.  This  error  can  be 
reduced  if  the  size  of  the  measurement  window  Is  enlarged  at  the  cost  of 
Increased  processing  time.  The  software  listing  for  this  study  is 
included  as  an  appendix. 


CONCLUSIONS 

1.  The  Hough  transform,  together  with  the  Sobel  edge  operator,  image¬ 
thresholding  methods,  and  decision  logics,  provides  an  effective  means  for 
detecting  and  recognizing  various  bridge  patterns  and  their  orientation 
angles  from  a  set  of  radar  imagery  of  the  Hudson  River,  New  York  State, 
area. 

2.  The  entire  set  of  the  selected  radar  bridge  image  patterns  was 
recognized  correctly  with  this  scheme  even  if  the  bridges  were  not 
centered  in  the  measurement  window. 

3.  The  resolution  of  bridge  orientation  angle  determination  with  a  32-  by 
32-pixels  measurement  window  is  approximately  22.5  degrees.  The 
resolution  can  be  increased  by  enlarging  the  measurement  window  size  as 
demonstrated  in  the  report. 

4.  The  quantization  error  of  the  measurement  window  can  also  be  reduced 
by  using  a  larger  window  size  at  the  cost  of  increasing  the  processing 
time. 


9 


e@ppee@@@@epppfi?pp@@ppppppe@pee@ppppppp@@p@pppppppppppp 

eepeppppppppppppepeepppppeeepeeeetieepeeeeppeppeeeeeeep 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 
pppppeppppeppeppppppppepppepeeeeeefipeeeeeepeeeeeeeeeee 
e@p@eppee@pe@pp@pppe@ppp@@pppppp@ppp@pp@p@pp@ppppppppp 
ppepeeppe@ppppepppe§eeppeppp@@ep@ppe@@e@eppp@@e@epeepp 
pppp  ppppppppppppppppppppppppppppppppppppppppppppppppp 

PP  PPPPPPPPPPPPPP  PPP  PPPPPPPPPPPPPPP  PPPPPPPPP 

pppppppppppp  pp  pppppppppp  pp  gp 

pp  ppppp 


gppppgpggppgppppppg 

pppppppppppppppppppppppp  ppppppppppppppppppp  @p 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 
ppppgppppppppppppppppppgpgppppgppppppppppgpgpppppppppp 
pppppppppeppppppppppppppppppeppepppppppeepppppeppeepep 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppppppppppeepppppppppppeppppeppeeeeep 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

ppppppppppeepppppppppepppppeeepppppppeppeppepepppeeeee 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

ppppppppppppppppppppppppppepepeeppepeppeppppeppepeeeep 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppepppppppppppppeepppepeepppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 


(a)  Pictorial  Print  of  Processed  Input  Image. 


I ROW  =  30 
ICOL  =  3 

AiNS  =  1 

ITT  -  22. 

J  =  32 

Ri  =  0.0 

Bridge  Detected  at  0.0  Degrees  with  X-Axis. 

FIGURE  2.  (b)  Recognition  Result  of  Bridge  Pattern  at  0  Degree. 
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pppppppppppppppppppppppppppppppppppppppppppppppppppppp 
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pppppppppppppppppppppppppppppppppppppppppppppppppppppp 
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(a)  Pictorial  Print  of  Processed  Input  Inage. 


IROW  =  9 
ICOL  =  7 
ANS  =  1 

ITT  -  22. 

J  =  31 

R1  «  22.5 

Bridge  Detected  at 

FIGURE  3.  (b)  Recognition  Result 


22.5  Degrees  with  X-Axis. 

of  Bridge  Pattern  at  22.5  Degrees. 
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pppppppppppppppppppppppp  ppppppppppppppppppppppppp 
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(a)  Pictorial  Print  of  Processed  Input  Inage. 


IROW  -  14 
ICOL  -  12 
ANS  *  1 

ITT  -  22 . 

J  -  30 

R1  -  45.0 

Bridge  Detected  at  45.0  Degrees  with  X-Axis. 

FIGURE  4.  (b)  Recognition  Result  of  Bridge  Pattern  at  45  Degrees. 
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(a)  Pictorial  Print  of  Processed  Input  Iaage. 


I  ROW  - 
ICOL  - 
ANS  = 
ITT  - 

J  - 

R1  * 
Bridge 
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30 

1 

22. 

29 

67.5 

Detected  at  67.5  Degrees  with  X-Axis. 


FIGURE  5.  (b)  Recognition  Result  of  Bridge  Pattern  at  67.5  Degrees 
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IROW  =  2 

ICOL  »  13 
ANS  =  1 

ITT  -  22 . 

J  =  28 

R1  =  90.0 

Bridge  Detected  at  90.0  Degrees  with  X-Axis. 


FIGURE  6.  (b)  Recognition  Result  of  Bridge  Pattern  at  90  Degrees 
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IROW  =  15 
ICOL  =  30 
ANS  -  1 

ITT  -  29. 

J  -  27 

R1  -  12.5 

Bridge  Detected  at  112.5  Degrees  with  X-Axis. 

FIGURE  7.  (b)  Recognition  Result  of  Bridge  Pattern  at  112.5  Degrees. 
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(a)  Pictorial  Print  of  Processed  Input  Inage. 


IROW  »  14 

ICOL  -  15 

ANS  -  1 

ITT  -  Id. 

J  -  26 

R1  -  135.0 

Bridge  Detected  at  135.0  Degrees  with  X-Axis. 

FIGURE  8.  (b)  Recognition  Result  of  Bridge  Pattern  at  135  Degrees. 
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FIGURE  8.  (c)  Hough  Transform  of  Bridge  Pattern  at  135  Degrees. 


23 


-^  -m  imTiiiraiftaa'ni  Hi  i  r  urniiiaBniiM  'aim  liiHrnHittr 


@eee@@@@@ee@e@@@@@@@eeeeee©@e@ee@@e@ee@e@eeee@ee@@@@e@ 

eeeeeeeeeeeeeeeeeeeeeeeeeeppppppppppppppppppppppeppppp 

eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee@eeeee@eeeeee©ee 

pepepeeeppppppppeepppppppppppppepppppppppppppppeeppppp 

ppppppppppppppppppppepppppepppeeeeppeeeeeeeeeeeepppppp 

ppppppeepeepppppeeppppppppppppppeppppppppppppppppppppp 

eeeeeeeeeeeeep@eeeeeeeeepeeee©©ee©©@ee@eeeee©©e©e©eeee 

pppppppppppppppppppeppppppppppppppppppppppppppppppepee 

eeeeeeeeeeeeeeeeeepeeeeeee@©@e@ee@eeeeeeeee@eeeeee©eee 

pppppppepppgeppppppppppppppppeppeppppppppppepppeeepppp 


peeeepeeeepppppppppppepppeppppppppeppppppppppppppppp 

eee@ggge@egee@pppgpp@i!>ees»@@@g@§@peepp@p§ppppt!> 
gg  ggggggggggggggeggeggpggggggggggpggpgggg 

gggp  gpepgggpgppppgppgpggppgppppgppppppe 

gggpgpggg  ppp  ggpggppegpeeppppggpgggpgp 

ggggggggggpe  ppppppppppppppppppp 

eggggpgpggppgpp  ppggpgppppgpggpppgg 

pgpggpggggpggpppggp  ppppppppppppppppppp 

gggggggggggggggggggggp  00  ggggggggg 

ggggggpggggggepgppppgppgp  ppp© 

gggpggggggpggpgpgpggggppppp  ©e 

ggpppggpggggpggppppppgpppppppgpppp  pp 

ggppgpgpggggggpggppgppppggppgpgggpggggp 

pppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppeepppppppppppppppeppppepppepp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

gggpgggggppppgpgepggpppgppgpggppgppppgppggpgpgggpgggpg 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

ppppepeeppepeppeepepppppeepeppepeepppppppppppppeepeepp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 


(a)  Pictorial  Print  of  Processed  Input  Image. 


IROW  =■  29 

ICOL  -  14 

ANS  »  1 

ITT  -  14. 

J  =*  25 

R1  =  157.5 

Bridge  Detected  at  157.5  Degrees  with  X-Axis. 


PIGORE  9.  (b)  Recognition  Result  of  Bridge  Pattern  at  157.5  Degrees 


pppeppppppppppppppeepeppppppppppppeppppppeeepppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppeeppppppppppppppppppppppppppppppppppppppppppp 

pppepppppppppppppppppppppppppppppppppppppppppppppppppp 

ppppppeesepppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

PPPPP9PPP8PPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPPP 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppppppppeeeeppppppppppppppppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

ppppppppppppppppppppppppppppppeepppppppppppppppppppppp 

ppppppppppeppppppppppeppppppeppppepppppppppppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 

ppeppppppppppppppppppppppppppppppppppppppppppppppppppp 

pppppppppppppppppppppppppppppppppppppppppppppppppppppp 


(a)  Pictorial  Print  of  Processed  Input  Inage. 


I  ROW  = 

0 

ICOL  = 

0 

ANS  = 

0 

ITT  - 

47. 

J  = 

25 

R1  =■ 

157.5 

"not  a 

bridge 

FIGURE  10.  (b)  Recognition  Result  of  an  Entirely  Black  Area. 


26 


333 


si xv  :<  h i  rn  s-nmnn  nir-iv’.  iv  (Hiniiw  ■»;»>! i rt” 


IT 


xxxaixxxaxxxxxxxxxa  xxxx  xxx®  xxx  xxx  xxxxxx  xxxxxx 

xxxxx® xxx®  igiiii  xxxxxxxxxxxxxxxxxxxx xxxx xxxxxx 

M’U.iiiiiittnii'itiitiiiitiitmittitiiiixiiiiitM 

MistiiiiiUiti.iiitttiitiiait'iitmiittiiiiii'i'i 

x-x  xxx  xxxxx  xxxx  xxxx  xxxx  xxxx  xxxx  xxxx  xxxx  xxxxx 

UMtiintutiiil'ittmitttmtxtiitutxttut 

M'Xt»tMX|>MitXltX»XXlXXaMlXt>XtX»X»XXXktXX«X 

Mm»*MtXXXXXl»XXXtXXXXXXXtttX*XiX*kXSXX»X 

xtmixittiiiiixxxxixiiitmiiimttxiixxx 

M;l*XtM«tttttXX*XXXlXXXaXX»tXl;tXtXXXXXXXX 

MnxxixtxixxtxxtxstxxMttxsiixstxtxxtxxtx 

mXlXMiXtXXXaXXXXX'tXXXtXXXXlXXXXXXXXtXXXt 

XXllKlIlHtStM'iXlXSttXXXM'XtXtXl'XXXXXt'XXlXX 

lAMltl'UlXXjiXXalXtXXl'XXlXtXttlXlXliXXXlXlX 

Xl'MlXJAiniXXXXtXt'XlXl’t'ttXXXXlXtXXXXU'XXl 

xxtxitutmmtxxixitixtxititmxxxxitxx 

XXXXX1XX1XXX1X1XXXXXXX11XXXX1XXXX1XXXXXXX 

XXXXltXXXXXHXXXXXXXXXXXXXXXXlXXXXXXXXXXX 

MXXXXXXXXXXXXXXXXtXXXXXXXXX'XXXXXXXXXXXX 

XlMXlXXXtM'Xl'lXXXXXXXXXXXXXXlXXXXXXXXXX 

XXXXXXXXXXXXXXXXXXXXXX1XX1XXX1XXXXXXXXX 

XXXXXXXXXXXXXIIXXXIXXXXXIXXXIXXXXXXXXXI 

XX1XXXXXXXXX1XXX1XXXXX1XXXX11X1XXXXXXXX 

ItXXXXXXXXXXXXlXXXXXlXXXXXXXXXlXlXXXXtl 

xxxxxxxxxtxxxxixixxxxxxixxxiixxixxxxxxx 

txxxxxxitxmxxxxxxxxxxxxxixxxxxxxxxxx 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

uxxxtxxxxxxxxiixxxxixxxxxxxitxxxxxixx 

XXXXXXXXXXXXXXXXIXXXIXXXXXIXXXXXIXXXX 

ixtxixxxxxxxixxxxxxxxxxxxxxxxxxxxxxxx 

XXXXXXXXXXXlXXXXXXXXXXXXXXX’XXXXXXXXXX 

xxxaxxxxxxxxxxxtxxxmxxxxxxxmxxxix 

IXXXXXXXXXXXXIIXIXXXIXIXXXXXXXXXXXXX 

XXXXXXMlSXXXXXXXXXlXXlXXXXXXXXXXllXX 

xxxxxxxxxtxxxxxxxxxxxxxxxxxxxxxxxxxx 

xxxxixuxxixxxxxxxxxuxxixxxxxxxtxxx 

XXXXIXIXXXX’XXXIXXXXXXXIXXXXXXXIXXIXX 

XXXXlXMiXXXXXXXXXXXXXXXXXXXlXXlXXXX 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXIXXIXXI 

xxxxxxxxsxaxxixxxxxxxxxxxiixxxxxtxx 

XXXXXXXXXXXXXXXXXXXXXX1XXXXX1XXXXX 
XXXXX'XXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
1XXXX  XXXXXXXX'XIXXXXXXXXXXXXXX  X  XXX'  X 
XXXXXXXXXXXXXXIXXXXXIXXXXXXXXXXXXX 
xxaxxxxxtxxxxiixsiixtxxxxxxiixixix 
xxxxxxxxxxxxxtixxxxxxxxxxzxxxxxxx 

1XXXXXXXXKXXXX1XXXXX1X1XXXXX1XXXX 

xtxxxixxtxxixxixxxxxxxxxxxxxxxixx 

xxaxsxxxixaxxxxxxxxxxxsxxtxxxxxx 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXIXXX 

XXXXlXXXXXXXXXXl'XlXXXXXXXXltXllX 

XXXXXXXXXXXXXXXXXXIXXX'XXIXXXXXIX 

XXXXIXXX'XXXIXI'XXXXXXIIXXXXXXXXXX 

ixxiixxxxxxxxxxxxxxxxxxxxxmxi 

axxxxxxxaxi'ixxxtxxxxxxitxxixxxi 

XXXXXXXXXXlXXXXltXlXlXXXlXXXlll 

xxxxx xxxxxxxxx xxxxxxxxxxxxxxxxx 

XltXXXXSXXXXXXXXXXXXXXlXXlXXXlX 
XXXXXXXUXXXXXXXXXXX1X1XXXXXX1 
XXXXIXXXXXXXXXXXI'XXXXXXXXXIXXX 
xxx a. x  xxxxxxxxa  xxxxxxx xx-xxxxax  x 
XXXXXXXIlXllXXXlXXlXXXl'iXXlXlXa 

xxxxxxxxixxixxxtxxtxixxxxxxxi 
xxxxxtxxxixtxxxxxxixxxxtxxxxi 
XXXXXXXXXXXXXXXXXlllXXXXXXXtX 
XXXXll'XXlXXXXXXXlXXXXXXXXXlXX 
XXX XX  XX XX XXXXXXXX X'XXIXXX'X XX X 
xxixxxxxxtxixxxxxixxxxxix'taxx 
XXXXXi  11XXXXXXXXXXXXX1X1XXXXX 
XXiXiXX  XXMXXXXXXXXlXl'XXXltlXXX 

istxix  xaxaixxxxxxxxxxxxxxxxxx 

XXXXXXXXXXXXXXXXXX'X  XXXXXXXX 

x  xx  xx  xxxxxx  xx  x-  x  xx  XX  a  x  a-  x  X  X  ax 
xxx XX  xxxxxxxxxxtxxxxxxixxx 
xxxixx1  auaixx  x  x  xx  x  x  xxx  xxx  i  xxx  x 

XXXXX  XXXX  XXXX  X  XXX XXXX X  X XXX 
XXXXX  XX xxxxxxx xxxxxx xxxxxx 
XXX  XXXXXX  XXXXXX  XXXXXXXX  xxx 
XXXXX  XXXX  XXXXX XXXXXXX XX XX 
XXXXX  xxxxxxxxxxxxxxxxxxxx 
XXXXX  XXXXXXXXXlXXXXXlXXll 

XXXXX  xxxxxxxxxxxxxxxxxxxx 

XXXXX  XXXX  XXXX  XX  xxxxxx  xxx 
XXXXX  XXXXXX  X  X  XXXXX XX XXXX 
XXX  XX  XXXX  XXXX  XXXX  XXX  XXXX 
XXXXX  XXXX  XXXXX  xxxxxxx  xxx 
XXX XX  XX XXXXXXX  XXX XX XX xxx 
XXXXX  XXXXXXXXXX'XXXXXXXX 
XXXXX  XXXXXXXXX’XXXXXXXXX 

XXXXX  XXXX  XX  xxxxxx  xxxxxx 
XXXXX  xxxxxxxxxxxxxxxxxx 
XXXXX  MjXXXXXXXXXXXXXXX'X 


XXXXX  xxxxxxxxxxxxxxxxx 
XXXXX  XXXX  XX XX XXXX  XXXX X 
XXXXX  xxxxxxxxxxxxxxxxx 
xxxxx  xxxxxxxaxxxxxxxx 
xxxxx  axxxxxxxxxxxxxxx 
xxx  xx  ax  xx  xx  xx  xx  xx  xx  xx 
xxxxx  axaxxxxxxxxxxxxx 
xxxxx  axaxxxxxxxxxxxxx 
JMia  aa.ax xxxx  xxxx  xxx 
xxx  x*  ax  ax  xx  xx  xx  xx  xxx 

XXXXX  XXXXXX XX XX XX xxx 

x  xx  xx  a X'  xx  xx  xx  x  x  xx  x  x 

»x« XX  XXXX  XX XX  XX  XX  XX 

a  am  xa  xx  xx  xx  xx  xx  am  x  x 
xxxxx  xx  xx  xx  xx  xx  am  xx 


X 
X 
X 
X 
XX 
XX 
XX 
XX 
XX 
X  XX 
xxx 
xxx 

XX  XX 
XXXX 
XXXX 
XXXX 
XXXX 

xxxxx 
xxxxx 
xxxxx 
xxxxxx 
xxxxxx 
xxxxxx 
xxxxxx 
xxxxxxx 
a  xxx! xx 
xxxxxxx 

XXXXXXX' 

xxxxxxx 

XXXXXXXX 

XXXXXXXX 

XXXXXXXX 

axxxxxxxx 
xxx xxxxxx 
xxx xxxxxx 

X XXXXXXXX 
X XXXXXXXX 

xxxxxxxx ax 
xxxaxxxxxx 
XXXXXXXX  XX 
XXXXX XXXX  XX 
xxx xxxxxxxx 
xuxxmin 

XXXXX XXXX XX 
XXXXXXXX  XX  X 
XXXXXX  XXXX 1 X 
XXXX xxxxxxxx 
XXXXXXXXXX  XX 

xxxx'xxxxxaxa 
xxxxxxxxxxxxx 
xxx  xx  xxxxxx  la¬ 
in  xxxxxxxx  ax 

2'XXX  XXXXXXX  XX 
XXXX  X'XX'XX  X  xxx 
xxxxxxxxxxxxx a 
xxxxx xxxxxx U 1 

XXXXXXXXXXXXX X 
XXXXXXXXXXXXX XX 
X  XXXXXXXXXXXXX  X 
XXXXXXXXXXXXX  XX 
XXXX xxxxxx xxxxx 
XXXX XX xxxxx 2 xxx 
XXXXXXX XXXXXXXX X 
XXXX1XXXXXXXXXXX 
X  XXXX  XX XXXXXXXX 1 
xxxxxxxxxxxxxxxxx 
XX xxx XXXXXX XXXX XX 
XXXXXXXXXXXXXXXXX 
XXXXXXXXXXXXXXXXX 

xxxxxxxxxxxxxxxxx 
xxx  x  xxxxx-  XXXXXX  XX  X 
xxxxxxxxxxxxxxxxxx 
xxxx xxx xxxxxxxx xx a 
xxxxxxxxxxxxxxxxxx 
xxxxxxxxxxxxxxxxxx  X 
XXXX XXX XXX XXXXXXXX 1 
XXXXXXXX  XX xxxxxx XX  X 
XX XX XXXX XX xxxxxx xxx 
xxxx xxxx XX xxxxxx xxx 
xxxxxxxxxxxxxxxxxxxx 
XXXXX  XX XX  XX  xxxxxx xx-x 
xxxxxxxxxxxxxxxxxxxx 
1 XXXXXXXX  XX  XXXXXXXX 1 
XXXXX XX XX XX XXXXXXX IX 
XXXXXXXX  XX  X  X xx xxxxxxx 
XX XX XX XXXX XX  XX xxxxxxx 
XXXXXX  XXXX  XX  XX xxxx xxx 
XXXXXXXXXXXXX  xxxxxxxxx 
XXXXX  XX  XX  XX  XX  XX  xx  xx  xx  x 
x xxxxxx xxxaxx xxxxxxxxx 
xxx  xxxx  xx  am  xx  xxxxxxxxx 
xxxxxxx xxxx  XX  xxxxxx xxx 
xx  xx  xx  xx  xx  am  xa  xxxxxx  xx  x 

XXXXXXXXXXXXX'X  XX  XX  XX  XX"  X 
XX  XX  XX  XX:  XX  XX  XX  XX  XX  XX  X  X  X 

x  xx  xxxx  xx  xx  xx  xx  ax  xxxx  xx  x 
xxxxxxx  xxxxxx  xa  xxxxxxxxx 
xxxxxxxxx  xxxx  xx  xa  xx  xx  xx  x 
x  xx  xx  xx  xx  xx  xx  xx  aa  xxxx  xx  x 


28 


Pictorial  Print  of  Processed  Input  Image,  and  Recognition  Result  of  Bridge 
Pattern  at  162  Degrees  with  Large  Measurement  Window. 


FIGURE  12.  Pictorial  Print  of  Processed  Input  Image ,  and  Recognition  Result  of  Bridge 
Pattern  at  172  Degrees  with  Large  Measureaent  Window. 
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C**«*****«****«**PROGRAM  -RD5- - REV  0S/14/H2********* *»»*»»****** 

C  IMIS  PROGRAM  RECOGNIZES  A  BRIDGE  OVER  WATER  USING  THE  HOUGH 

C  IRANSFURh  WHICH  CONVERTS  A  LINE  IN10  A  POINT.  THE  EDGES  ARE 

C  DEI  EC  I  ED  USING  A  80BEL  OPERATION  (SUBROUTINE:  *SOBEL“>. 

C  MODIFIED  TO  RUN  IN  THE  RTE  OP.  SYS . 

C 

C 

C 

C 

PROGRAM  RD5 

COMMON  IXPOS . I YPOS , IRUN , I MAGE . IAi 
COMMON  LEXLU 
DIMENSION  IA3< 1026) 

EWIII  VALENCE  < 1 A1 . IA3<  3 ) > 

DIMENSION  IA1<1024) , I A2( 1 U24 ) 

DIMENSION  IEDGE ( 1024 ), IDCB < 144) , IXBUF1  4?) , ITBUF ( 4?  > 

DIMENSION  IXFRM(64 ,32) , IA<64 ,32 > ,MAX< 1 1 > , I  MAXI  11 ) 

DIMENSION  JMAX(ll) 

DIMENSION  I  DARK (512  > ,INPUT(43> .LABIN(iO) .1  ABPT(IO) ,L.APEG( 1 0 ) 
DIMENSION  NAME (3) .LABHIH  10) ,LABA( 10) 

DIMENSION  IWKD(7) ,1WR2(7) 

DA  I A  INPOT/2HIS.2HIN.2HPS.2HPI .2HJP.2HTR.2HSR .2HI1 , 2HE V , 2H 1 H , 2HF  D , 
1 2H0R , 2HLU . 2HX V . 2HRS . 2HNI . 2HDK , 2MNA , 2HIM , 2HRU . 2HR A . 2HR I . 2HCR . 2HP A , 
22MGL . 2HMA , 2HMN , 2H  TF . 2HST , 2HHG . ?HHS , 2 HAP , 2HSB . 2H0E , 2HPE , 2HHU , 2HPH , 
32HT  T . 2HBD . 2HTS , 2HLX . 2HFA . 2HF  A/ 

DATA  LABIN/2HIN , 2HPU ,2HT  ,7*2H  / 

DATA  LABPT/2HIN.2HPU.2HT  , 2HSC , 2HAL , SHED , 4*2H  / 

DATA  I YES/2HYE/ 

DATA  LABEG/2HED ,2HGE  ,80211  / 

C 

C  SET  VARIABLES  TO  DELAULT  VALUES 
C 

IC=0 

CALL  RHPAR(IAl) 

IF  (  IAI  .  EO ,  0  )  I  Al=*l 

Ltl-IAl 

LUOT-LU 

CALL  LUNT (75B .LUWAT . I > 

IT  1*0 

3  WR I TE  (  LIT  1 2 ) 

2  FORMAT ( "STEP  FILE  NAME?") 

RE AD (LU, 30) NAME 

4  CALL  OPEN! IDCB.I. NAME) 

CALL  FEROR( I ,4 ,LU> 

IF  (T.LT.O)GO  TO  3 

CALL  READF ( IDCB . IERR . IXBUE ,47 , IMAC9) 

CALL  READF ( IDCB , IERR , IYBUF , 4? , IMAG3 ) 

CALL  CLOSE (IDCB) 

IMAGS-IMACS+1 
WRITE (LU . 1 ) IMAG8 

1  FORMA T( "NO.  OF  IMAGES-", 12) 

NAV-S 
IXPOS-O 
1 VP 08=0 
IMAGE=t 
J  RUN-0 
NCR AY-16 
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pr 


oosv 

1  111*0 

ni)60 

ISS*1 

0061 

IDIR=1 

006? 

I  OR*? 

0063 

111=0 

0064 

HAX  =  0 

0065 

MIN=0 

0066 

NCHRS-7? 

0067 

1DIM*32 

0068 

IA3(  1 )=HAX 

006V 

I  A3 ( 2 ) =HIN 

0070 

IB*-1 

0071 

C 

007? 

C 

INl’Ur  COMMAND  LOOP 

0073 

C 

0074 

C 

DK  — INPUT  ARRAY  DARK  LEVELS 

0075 

c 

LU— INPUT  LIST  l.U 

0  076 

c 

IN— INPUT  ARRAY  DATA 

0077 

c 

IT  — INPUT  THRESHOLD  CONSTANT  TOR  INPUT 

IHAGE 

0078 

e 

PI— PRINT  THE  INPUT  ARRAY 

0079 

c 

TR— TERMINATE  PROGRAM 

0080 

c 

TH— THRESHOLD  FOR  BACKGROUND 

0081 

c 

SP — StIBEL  OPERATION 

008? 

c 

OE— OUTPUT  RESULT  OF  SOBEL  OPERATION 

0083 

c 

HU— HOUGH  TRANSFORM 

0084 

c 

PH— OUTPUT  HOUGH  TRANSFORM  MATRIX 

0085 

c 

TT— THRESHOLD  VALUE  FOR  DETECTING  AN  EDGE 

0086 

c 

BD— BRIDGE  DETECTION 

0087 

c 

TS— INPUT  THRESHOLD  VALUE  BEFORE  SOBEL 

OPERATION 

0088 

c 

LX— DISPLAY  RESULT  OF  BRIDGE  DETECTION 

ROUTINE  ON 

0089 

c 

PE— DISPLAY  HARD  COPY  OF  INPUT  IMAGE 

0090 

c 

XY— SLEW  TO  AN  IMAGE 

0091 

c 

IS— INPUT  STAGE  STEPS  FOR  NEXT  IMAGE 

0092 

c 

OR  — INPUT  ORIENTATION  VALUE 

0093 

c 

RS— RESET  TO  FIRST  IMAGE 

0094 

c 

NI— SLEW  TO  NEXT  IMAGE  IN  BEOUENCE 

0095 

c 

NA— INPUT  NUMBER  OF  FRAMES  AVERAGED 

0096 

c 

IM— SLEW  TO  SELECTED  IMAGE 

0097 

c 

RU— AUTOMATIC  RUN  THROUGH  IMAGES 

0098 

c 

RA— RUN  AGAIN  STARTING  AT  STATEMENT  3 

0099 

c 

R I— REMOVE  6ELECTED  IMAGE 

0100 

c 

0101 

10 

WRITE1LU.20) 

010? 

20 

FORMAT!*??") 

0103 

READ1LU .30 ) ICMND 

0104 

30 

FORMAT (3A2> 

0105 

IF  (ICMND. EQ. INPUT! 1)>G0  TO  100 

0106 

IF  ! ICMND. EQ.INPUT(2)>GO  TO  200 

0107 

IF  <  ICMND. EO.  INPUT  <B>  >READ(t.U,«) IT  1 

0108 

IF  (ICMND. E0.INPUT!4))G0  TO  300 

0109 

IF  ( ICMND. EQ.  INPUT!  10  >  >READ(LU ,» >'ITH 

0110 

IF  (ICMND. EQ.INPUr(ll))GU  TO  6000 

Dili 

IF  (ICMND. EQ.INPUT(20))GO  TO  1300 

on? 

IF  ( ICMND. EQ. INPUT! 15) ICALL  RESET 

0113 

IF  ( ICMND. EQ. INPUT < 16) >G0  TO  1000 

0114 

IF  ( ICMND. EQ. INPUT! IS) >READ(LU.*)NAV 

0115 

IF  ! ICHND . EQ . INPUT! 19) )GO  TO  1200 

0116 

IF  ! ICMND. FQ.INPUT(21))GO  TO  2100 

0117 

IF  (ICMND. EQ.INPUT!2?))GO  TO  2300 

0118 

IF  (ICMND. EQ. INPUT! 17) >GO  TO  1100 
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Oil? 

IF  ( ICMND . EQ . INPUT ( 13) >GO  TO  500 

01. ?0 

IF  < ICMND. Ed. INPUT<6))GO  TO  1400 

0121 

IF  < ICMND . EQ. INPU T(33) )  GOTO  4000 

0122 

IF  < ICMND . EQ. INPUT <3S> )  GOTO  4200 

0123 

IF  < ICMND .EQ. INPUT! 371)  COTO  4600 

0124 

IF  < ICMND. EQ. INPUT! 34 >>  GOTO  4100 

012S 

IF  < ICMND . EQ. INPUT <361 1  GOTO  4500 

0126 

IF  <  ICHND .EQ.  INPUT <381  >READ(l.U ,  *  >  ITT 

0127 

IF  < ICMND .EQ. INPUT <141 (GOTO  600 

0128 

IF  < ICMND. EQ. INPUT! 40) )READ<LU, *> IT2 

0129 

IF  < ICMND. EQ. INPUT! 41)1  GOTO  987 

0130 

IF  ( ICHND . FQ. INPUT <241 >  GO  TO  6800 

0131 

31 

FORMAT  1 *1*1 

0132 

CO  TO  10 

0133 

C 

0134 

C 

STORE  STAGE  STEPS  FOR  NEXT  IMAGE  OR  INSERTED 

IMAGE 

0135 

C 

0136 

C 

SLEW  TO  STARTING  IMAGE.  0  INPUT  ABORTS 

0137 

c 

SECTION. 

0138 

c 

013? 

c 

0140 

100 

UR  I  TE(LIJ  .101) 

0141 

101 

F OR MAT < “START  AT  IMAGE  0?"> 

0142 

READ1LU,*)! 

0143 

IF  <  I .  GT  . I MAGS .OR . I . LT . 1 )GO  TO  10 

0144 

IRUN-2 

0145 

CALL  STEP  S  <1 XBUF . I YMIF , I MAGE , I , 1XSTP , IYBTP  > 

0146 

IXPOB-IXPOSTIXSTP 

0147 

1  YPOS-IYPOB-UYSIP 

0148 

I MAGE- I 

014? 

CALL  SSLEM(-IXSrP ,-IVSTP , 1 ) 

0150 

c 

0151 

c 

INPUT  ADD  OR  INSERTIMAGE  FLAG.  "ICMND". 

AND 

HOVE  TO 

0152 

c 

NEU  IMAGE  WITH  KURD  CONTROL. 

0153 

c 

0154 

WRITE1LU ,120) 

0155 

120 

FORMAT < "ADD( 1 )  OR  INSERT<2)  IMAGES? *> 

0156 

READ<LU.«)ICMND 

01S7 

107 

WRITEILU. 1021 

0158 

102 

FORMAT ( "MOVE  STAGES  WITH  KBRD  UNTIL  DESIRED 

IHACE  IS") 

0159 

WRITEILU, 103) 

0160 

103 

FORMAT < "DISPLAYED .  ENTER  1  TO  STORE  NEW 

IMAGE  POSITION") 

0161 

WRITEILU. 109) 

0162 

10? 

FORMAT < “ENTER  2  FOR  FEATURE  DETECTION") 

0163 

IXD-0 

0164 

IYD-0 

0165 

105 

RE AD  <  LU . * ) I XSTP . I VSTP 

0166 

IFdXSTP.EO.  1>G0  TO  104 

0167 

IF <IXSTP.EO.2IGO  TO  200 

0168 

CALL  SSLEWl -1X8  TP , -IYSTP , 1 ) 

0169 

IXD=IXD+IXSTP 

0170 

I YD-I YD+IY8TP 

0171 

GO  TO  105 

0172 

C 

0173 

C 

COMPUTE  NFW  STAGE  POSITIONS  RELATIVE  TO 

STARTINC 

0174 

C 

IMAGE  AND  BRANCH  TO  ADD  INSERT  SECTION. 

0175 

C 

0176 

104 

txras-ixPos+ixD 

0177 

IYPOS-1 YPOS+IYD 

0178 

It  (  ICMND. FO.?)GI)  TO  1  1  1 

3H 


ADD  IMAGE.  SECT  ION  i  STORE  NEXT  IMAGE  STEPS 


017? 

0180 

0181 

0102 

018J 

0184 

0185 

0186 

0187 

0188 

018? 

01?0 

01?1 

01?2 

0193 

0194 

0195 

0196 

0197 

0198 

0199 

0200 

0201 

0202 

0203 

0204 

0205 

0206 

0207 

0208 

0209 

0210 

0211 

0212 

0213 

0214 

0215 

0216 

0217 

0218 

0219 

0220 

0221 

0222 

0223 

0224 

0225 

0226 

0227 

0228 

022? 

0230 

0231 

0232 

0233 

0234 

0235 

0236 

0237 

0238 


C 
C 
C 

1XBUFI  IMAGE  >«IXD 
lYBOF'C  I NAGE)-I YD 
GO  10  114 
C 

C  INSERT  IMAGE  SECTION!  INSERT  NEU  STEPS  INTO  CURRENT 

C  IMAGE  INDEX,  MOVE  I  XBUF  AND  IYBUF  DOWN.  COMPOTE  STEPS 

C  FROM  NEU  IMAGE  TO  OLD  NEXT  IMAGE  AND  STORE  IN  INDEX41 

C 

111  1 XS  TP -I XBUF  < IMAGE  >-IXD 
IYS1P-IYBIIF(IMAGE>-IYI> 

I XBUF  < IMAGE  >  =  IXD 
IYBUF < IMAGE  >-IYD 

IF (IMAGS- IMAGE. LT. 2 >GO  TO  115 
DO  113  1-1 , IMAGS-IMAGE-1 
I  XBUF ( I MAGS- I ♦ 1 > - I XBUF < I MAGS- I > 

113  IYBUF  < IMAGS- 1 4l >-IYBUF  < I MAGS- 1 > 

115  IXBUF1 IMAGE+l )-IX8TP 

IYBUF  < IHAGE+1 >*IYSTP 
IMAGS-1MAG841 
C 

C  ASK  IF  10  CONTINUE  AND  IF  NOT  STORE  STEP  FILE 
C  IN  DESIGNATED  NAME 

C 

114  1 MAGE ~ I MAGE *  1 

IF ( IMAGE . EO . 50 .OR . IMAGS .EQ. 50 >GU  TO  108 
UR  ITE  (1.(1 , 1 06  > 

106  FORMA T( "CONTINUE? *> 

READ(l.U  ,30)1 
IF  < I .EQ . IYE8)G0  TO  107 
108  IF(ICMND.E0.1)IMAGS-IMAGE 

112  UR  ITE  (HI.  2) 

RE AD <LU, 30) NAME 
CALL  OPEN! IDCB, I .NAME) 

IF  ( I . EO. -6)CALL  CREAlIDCB, I, NAME, 1 .3.0,9) 

CALL  WRITF( IDCB , I , I XBUF , I MAGS- 1 ) 

CALL  UR I1F ( IDCB , I , IYBUF , I MAGS- 1 ) 

CALL  CLOSE! IDCB) 

CALL  RESET 
GO  TO  10 
C 

C  INPUT  AND  AVERAGE  ARRAY  FRAMES,  SUBTRACT  DARK 
C  LEVELS.  THRESHOLD  AND  ORIENT 

C 

200  CALL  AVGIN< IA1 , IA2 , i 024 ,NAV, I DARK ) 

C 

C  (HE  FOLLOWING  STATEMENTS  FIX  AN  APPARENT  HARDWARE  ERROR  BY  SORTING 

C  SUCH  THAT  IHE  FIRST  ROW  IN  *IA1"  BECOMES  THE  LAST  ROW. 

C 

DO  203  I»! .1024 
203  IA2(I)-IAt(I> 

DO  201  1=1 ,??2 

201  IA1<I>-IA1(  (432) 

DO  202  1-993,1024 

202  1 A1 ( I ) -IA2< 1-99? ) 

CALL  IOHNM  IA1 ,32, TOR) 

C  SE  f  THRESHOLD  FOR  INPUT  ARRAY 
O 


023? 

CALL  MATH  1 < IA1 , ITM) 

0340 

DO  220  I*=l . 1024 

0241 

IRIAKI)  ,(;E.ITH)XA1<I)-1 

0242 

220 

CONTINUE 

0243 

1MIRUN.GT.  0)4000. 10 

0244 

C 

0245 

C 

PRINT  INPUT  IMAGE 

0246 

C 

0247 

300 

CALL  P DATE (LUOT  > 

024B 

CALL  MTUUTdAl .LAPIN ,32 ,32 , 3 .LUOT , -NCHR3) 

024? 

GO  TO  10 

0250 

C 

0251 

C 

CHANGE  PARAMETERS 

0252 

C 

0253 

C 

INPUT  LIST  LU 

0254 

C 

0255 

500 

READ1LU , *  M.UOT 

0256 

IT  (LUOT  .  EQ .  6)NCHRS=*132 

0257 

IE  (l.UOT  .  EQ .  l.U)NCHRS=72 

02SB 

GO  TO  10 

025? 

C 

0260 

C 

INPUT  AND  PACK  DARK  LEVELS 

0261 

C 

0262 

1100 

CALL  AVGIN( IA1 , IA2, 1024 ,NAV> 

0263 

GO IU  10 

0264 

C 

0265 

1400 

STOP 

0266 

C 

0267 

C 

SLEW  TO  NEXT  IMAGE 

0266 

c: 

026? 

1000 

IE( IFBRK ( I > .LT . 0 )GO  TO  1010 

0270 

IT ( IMAGE. EU.IMAGSIGO  TO  1010 

0271 

CALL  6SLEW(-IXBUF (IMAGE) ,-IYPUF( IMAGE > , 1 > 

0272 

I XPOS« I XPUS* I XBUF ( IMAGE ) 

0273 

I YPOS* I YPOS4 1 TBUF ( IMAGE  > 

0274 

I MAGE -I MAGE *1 

0275 

IE(IRUN.EO. 1)200.10 

0276 

1010 

CALL  RESET 

0277 

GO  TO  10 

0278 

c 

027? 

c 

SLEW  TO  SELECTED  IMAGE 

0280 

c 

0281 

1200 

READ(LU.«)I 

0282 

I F ( I , G  T  . I MAGS . OR . I . L  T . 1 ) GO  TO  10 

0283 

CALL  STEPS (IXBUF , IYBUF , IMAGE , I , IXSTP . IYSTP  ) 

0284 

IXPOS=IXPU8+IX6TP 

0285 

I  YP08=*I  YPOSdYSTP 

0286 

IMAGE- I 

0287 

CALL  SSLEW(-IXSTP, -IYSTP, 1> 

0288 

GO  TO  10 

028? 

c 

02?0 

1300 

CALL  RESE1 

02?  1 

JRIINM 

02?2 

GO  TO  200 

02?3 

c 

02V4 

r. 

0295 

2.100 

URI IE(LU, 23101 

0296 

2310 

F  UR MAK  "IMAGE  0  TO  PE  REMOVFDT"  > 

02V/ 

R1  AD(l.U .  *  >  I 

02?M 

If (t.Ur.IHAGS.OR.I.LT.l)GO  TO  18 

36 


r 

■  I 


029V 

IF ! I . EQ. IHAGS )GO  10  2320 

0300 

CALL  STEPS! IXBUF , I YBUF , 1-1 , 1*1 , IXD , IYD) 

0301 

IXBUF !I-1)»IXD 

0303 

1YBUF !I-1)»IYD 

0303 

DO  2330  IXD»l,IHACS-2 

0304 

I XBUF  < IXD  >  =  IXBUF ! IXD-U > 

0305 

2330 

IYBUF1 IXD)*IVBUF! IXD*1 ) 

0306 

2320 

1 MAGS* IHAGS- 1 

0307 

MR  1 T  E ( LU , 2 ) 

0308 

READ1LU , 30 )NAME 

0309 

2350 

CALL  OPEN! IDCB , IXD .NAME ) 

0310 

IF  < IXD . EM . -6) CALL  GREAT < I DCS , IXD .NAME . 1 .3 . 0 .9) 

0311 

CALL  MRITF! IDCB , IXD , IXBUF , I ) 

0312 

CALL  MR ITF  < IDCB , IXD , IYBYF , I ) 

0313 

CALL  CLOSE! IDCB) 

0314 

GO  TO  10 

0315 

C 

SOBEL  OPERATION  i  USED  TO  DETECT  EDGES 

0316 

C 

0317 

4000 

CALL  6QHEL! IA1 , IEDGE ,32,32) 

0310 

C 

0319 

C 

SET  THRESHOLD  FOR  SOBEL  ARRAY  'IEDGE* 

0330 

C 

0321 

CALL  MATR2! IEDGE. IT2) 

0322 

IDIN*32 

0333 

IFORUN.GT.  0)4500, 10 

0324 

C 

032S 

C 

OUTPUT  RESULT  OF  SOBEL  OPERATION  TO  LINE  PRINTER 

0336 

c: 

0337 

4100 

CALL  MTOU 1 ! IEDGE . L ABEG . I DIN , I D I N , 3 , LOOT , NCHRS ) 

0328 

GOTO  10 

0329 

c 

0330 

c 

P ICTOR IALLY  PRINT  OUTPUT  OF  INPUT  ARRAY 

0331 

c 

0332 

4200 

CALL  LP IM! IA1 ,32) 

0333 

GOTO  10 

0334 

C 

0335 

C 

HOUGH  TRANSFORM  I  USED  TO  CONVERT  LINE  INIO  POINT 

0336 

t: 

0337 

4500 

1*1 

0330 

4501 

CALL  HOUGH !IEDGE,32.32,IT2,IXFRH.64.32> 

0339 

c 

0340 

IFIIRUN.GT. 0)6000, 10 

0341 

c 

0342 

c 

PRINT  HUUGM  TRANSFORM  MATRIX  ON  LINE  PRINTER 

0343 

u 

0344 

4600 

CALL  Ml OUT ! IXFRM , LABHU , 64 , 32 , 3 , LUOT , -NCHRS  > 

0345 

GO  ID  10 

0346 

C 

0347 

C 

SI.EM  TO  AN  IMAGE 

0340 

C 

0349 

600 

RF.AD(LU ,  0  >  IXS1P  .  I YSTP 

03S0 

IF ! IXSTP . EO. 1 >G0  TO  10 

0351 

CALL  S6LEMT- IXSTP , -I YSTP , 1 ) 

0352 

COIO  600 

0353 

987 

1C-IC+1 

0354 

1B»IB«1 

0355 

C 

0356 

C 

DISPIAY  THE  INPUT  ARRAY  'IAI*  ON  THE  LEXIDATA.  TH 

1 357 

c 

CAPABLE  IIF  DISPLAYING  12  ARRAYS,! 4  COLUMNS..)  ROMS 

0  Tr,U 

i: 

1IIKM  GEODE  Nil  ALLY  STARTING  AT  TUP  im  r  no  ill d  Am; 

DISPl  AY  IMF  INPUT  ARRAY  "IAI*  ON  THE  LEXIDATA.  THE  SUBROUTINE  IS 
CAP  AWE  IIK  DISPLAYING  12  ARRAYS, (4  COLUMNS  .  I  ROMS).  IT  DISPLAYS 
THEN  f,EIJUE  N I  1  ALLY  STARTING  AT  TOP  LEFT  CORNER  AND  ENDING  AI  THE 


0JS9 

0360 

0361 

0362 

0363 

0364 

•  36S 
0366 
0367 
0368 
0369 
0370 
0371 
0372 
037J 
0374 
0375 
0376 
0377 
0378 
0379 
0380 
0381 
0382 
0383 
0384 
0305 

•  386 
0387 
0388 
0389 
0390 
0391 
0J92 
0393 
0394 
0395 
0396 
0397 
0398 
0399 
0400 
0401 
0402 
0403 
0404 
0405 
0406 
0407 
0408 
040V 

•  410 
0411 

•  412 
0413 

•  414 
0415 
0416 
0417 
04 1 M 


C  HO non  RICH!  CORNER ,  l UK  WORD  “BRIDGE"  APPEARS  IF  A  BRIDGE  OVER 

C  WATER  IS  DETECTED.  At.SU ,  1HE  ANGLE  THE  BRIDGE  HAKES  WITH  THE  POB- 

C  ATIVE  X-AXIS  IS  SHOWN  ON  THE  LEXIDATA. 

C 

IF ( ANS . EQ . 0 )GO  TO  10 
I DEG-R 1 

CALL  LXKDT  IA1  . 1  DEC  ,  IC  ,  IB  ,  ANS> 

GO  TO  10 
C 

C  CLASSIFICATION  OF  BRIDGES 

C 

6000  CALL  BRDGE  < I A 1 , ANS , I COL , I ROW) 

UR I IETLUOT ,225) IRDM 

225  FORMAT! IX, “I ROW-  ■1I3> 

WRITE (LOOT .226) I COL 

226  FORMAT! IX, “ICOL=  *  ,  13) 

WRITETl.llOT  ,227)ANS 

227  FORMAT < IX , “ANS*  *.14) 

C  SET  THRESHOLD  FOR  HOUGH  TRANSFORM  ARRAY  “IXFRM* 

CALL  MATR3T IXFRM. ITT) 

WR ITETLUOT , 4S05) ITT 
4S0S  FORMAT (IX, "ITT*  ",F3.1i 
DU  60S0  J=1 ,32 
DO  6050  I « 1,64 
IA<I,J)*IXFRM<I.J> 

6050  CONTINUE 

DO  6700  K*1  ,10 

MX*0 

IMX=0 

JMX*0 

DO  6600  J-l ,32 
DO  6600  1=1 ,64 
1F(IA(I,J).LT. MX )GO  TO  6601 
MX»IA( I , J) 

IMX*=  I 
JMX*J 

6600  CONTINUE 
MAX<K)=MX 
TMAX1K  >  =  IMX 
JMAX(K)=JMX 
IA( IMX , JMX ) =0 
6700  CONTINUE 
6900  M=0 

DO  6910  K»1 ,10 

IF ( JMAX ( K ) . NC .JMAX < 1 > )GO  TO  6910 
lFTHAXTK)  .IF.  IT DGO  TO  6910 
H=M+1 

IF (M.LT ,2>G0  TO  6910 
J~JHAX(I ) 

GO  TO  6100 
6910  CONTINUE 

DO  6925  1*2,10 
M-0 

DO  6920  K-2,10 

I F <  JMAX  < K ) . NE . JMAX  < I ) ) GO  TO  6920 
IF  THAXTK  >  .  l.E  .  I  TT  >GO  TO  6920 
M-M*l 

IF1M.LT  ,2>r.n  TO  6920 
.1  =  JMAX!  I) 

GO  111  6100 


38 


0419 

6920 

tUNIINUF. 

0  420 

6925 

CONI  INIJE 

0421 

.1  -  JMAX  <  1  ) 

04.?;? 

no  it)  6100 

04?  3 

r 

0424 

6UU0 

CALI.  HI  lllll  <  I  A,  LARA,  64 ,  32, 3, LUO T  ,-NCHRS) 

0425 

6100 

UR  11 E (LUO! ,6150  >J 

0426 

6150 

f  ORMATdX,  “J*  ",IJ> 

0427 

R  1  - !  32- J  >  *22 . 5 

0428 

IF (J  .11 .20  >R 1“1 80  - J*2 . 5 

0429 

URI  TEILUOI  ,6160  > R 1 

0430 

6160 

F OR MAI ( 1 X , "R 1  -  ",F5.I ) 

0434 

IF  CANS.  EO.  1 . 01URI  IF.ILIIUI  ,6170  )Ri 

04.'  ? 

IF(ANS.rq.O.O)URI1MI.UOT.6ieO> 

0433 

6170 

FORMA! < IX. “PKIDGE  DETECTED  AT  “,F5.1.“  DECREES  Willi  X-AXIS") 

0434 

61811 

FORMAT! IX, "1HIS  19  NOT  A  PR  I  DOE  "  ) 

0  4,55 

) F  < IRUN . EQ  .  1 )GO  TO  6190 

0436 

lEdRUN.EQ.  2)105, 10 

0  437 

6190 

CALL  EXEC  <  1  ,  l.UWAT  ,150,1) 

0438 

GO  Tl)  1000 

0  439 

C 

RESET  THE  MAIN  PROGRAM  AND  INITIAL  VALUES 

0440 

t: 

0441 

2100 

CALL  RESET 

0442 

CO  Til  3 

0443 

INI) 

0444 

c 

0445 

<; 

IMIS  IS  A  SUPRUUTINE  10  RESET  THE  MAIN  PROGRAM  "RD“  TO  ITS 

0446 

c 

INITIAL  S 1  AGES 

0447 

c 

0448 

SlIPROUT  INE  RESET 

0449 

COMMON  IXP0S.IYP09, IRUN. IMAGE 

0450 

I  MAGE  =  l 

0451 

JRUN'O 

045? 

CALL  SBLEUIIXPOS.IYPOS.l) 

0453 

IXPOS=0 

0454 

I YPOS^O 

0455 

RETURN 

0456 

END 

0457 

c 

0458 

i: 

SUPRUUTINE  TO  CALCULATE  THE  NUMBER  OF  STEPS  BETWEEN 

0459 

c 

IMACES  IHI  AND  IM2  FROM  THE  SUCCESIVE  DISTANCES 

0460 

c 

STORED  IN  "IXBUF “  AND  "IYBUE* . 

0461 

c 

0462 

SUBROUTINE'  STEPS!  IXBUF  ,  IYBUE  ,  INI ,  IhZ  ,  IX8TP  ,  IYSTP  > 

0463 

DIMENSION  IXBUF !1 >, IYBUE! 1 > 

0464 

IXBTP-0 

0465 

IYSTP^O 

0466 

IF  ( INI  .  EIJ .  IM2)RETIJRN 

0467 

M-IM1 

0468 

N-IM2 

0  469 

IF! IMi .LT. IM2)GU  TO  20 

0470 

M»TM2 

0471 

N«=IMi 

0472 

?o 

DU  30  I*M,N-1 

0473 

IXSTP-IXPlIFd  >*IXSTP 

0474 

30 

IYSTP-IYBUF  (D  +  IYSTP 

0475 

IFUm.LT '.  IH2)RETURN 

0476 

TX31P  —  IX9IP 

047/ 

IYSTP— IYSTP 

0478 

PE  IIIRN 

39 


047V 

0400 

0401 

0403 

0403 

0484 

0485 

0486 

0407 

0408 

040V 

0490 

0491 

0492 

0493 

0494 

0  49S 

0496 

0497 

0  490 

0499 

0S00 

0501 

0S02 

0903 

0904 

0909 

0906 

0507 

0500 

0509 

0510 

0511 

0512 

0513 

0514 

0515 

0516 

0517 

0510 

0519 

0520 

0521 

0522 

0523 

0524 

0525 

0526 

0527 

0529 

0529 

0530 

0531 

0532 

0533 

0534 

0535 

0536 

0537 

ns  in 


END 

C 

C  AIM  OM  A 1 1C  THRESHOLD  FOR  INPUT  IMAGE.  THE  ARRAY  IS  CALLED  •IAf. 

C 

SUBROUTINE  MATR 1 < I A 1 . 1 T  H  > 

DIMENSION  IA1 < 1024) 

SUM«=0. 

1MAX-IA1 ( 1 > 

IXI  28  I -2. 1024 

IF  II Ail  I)  .  GT  .  I  MAXI )  IMAX1-IA1  If) 

28  CON I I NOE 

I TH*  .  500FI  OAT <  IHAX1 ) 

DO  38  1=1 .1024 
IF  < ITH. GT . 1AI < I ) )  I A1 ( I ) =0 
38  CONI INUE 
RE1URN 
END 
C 

C  SET  THRESHOLD  FOR  SOBEl.  OPERATION.  U3EB  ARRAY  "IEDGE* . 

c: 

SUBROUTINE  MATR21 IEDGE . IT2) 

DIMENSION  IFDGE! 1024) 

SUM=0 . 

IMAX=IEDUE < 1 > 

DU  28  1=2.1024 

IF  < IEDGE < I > .GT . I MAX) TMAX- IEDGE < I > 

28  CONTINUE 

M2=EL0AT(IMAX)*.SS 
DO  38  1=1,1024 

IE1IT2.GT . IEDGE!  I)  HEDGE!  I  >«0 
38  CONTINUE 
RETURN 
END 
C 

C  SET3  THRESHOLD  EOR  THE  HOUGH  TRANSFORM.  ARRAY  NAME t • IXFRM" . 

C 

SUBROUTINE  HATR31 IXFRM, ITT) 

DIMENSION  IXFRM! 64 . 32) 

IMAX=IXFRM! 1,1) 

Dll  20  I>1.64 
DO  28  J=1 , 32 

IE! IXFRM! I.JI.GT.IMAX) IMAX=IXFRM! I , J  > 

28  CONTINIIF 

ITT-IMAX0.75 

RETURN 

END 

SUBROUTINE  LXRD! IAI , IDEG , IC , IB .ANSI 


••••OOTHIS  DISPLAYS  A  32X32  ARRAY  AND  THE  RESULT  IIF  A  BRIDGE  DETECTION 

ooooooroutine 0000000000000000*0000000000000000000000000000000000000000 


COMMUN  LF.XI.U 

DIMENSION  IAI ! 1024) , INUM! 2 1 , 1WD! 31 
DATA  IUD/2HBM .2HID.2HCE/ 

1UNES-MUD! IDEG, 101 
IDEG=! IDEG-IIWESI/IO 
I IEN8=MUD! IDEG , 10  > 

IDEG=! I DEG- 1 TENS  >710 
I IIUM=MOD<  IDEG  ,10) 

I MUM! 1 >  »2S6*< IHUN» 4BM-! I1FN9M8> 


0539 
I)  540 
0541 
1)54? 
0543 
0544 
0545 
0546 
0547 
0548 
0549 
0550 
0551 
055? 
0553 
0554 
0555 
0556 
0557 
0S58 
0559 
0560 
0561 
0562 
0563 
0564 
0565 
0566 
0567 
0568 
0569 
0S70 
0S71 
0572 
0573 
0574 
0575 
0576 
0577 
0578 
0579 
0580 
05B1 
0582 
0583 
0584 
0585 
0586 
0587 
0588 
0589 
0590 
0591 
0S92 
0593 
0594 
0595 
0596 
0597 
0598 


C 

c 

c 


c 

c 

c 


25 


C 

C 

C 

C 

c 


10 


INUH<  2 ) -256*  < IONES+48 ) 

IR-IC-l 

il=o 

IR-IC-1 

JKIC.Gl  .4>IL=J 
IFdC . GT . 8) IL-2 
IF ( IR . EG . 4 ) IR=0 
IF  dR , EO . 5) IR* 1 
I F  < I R . EG . 6 ) IR-2 
IF(  IR  .  EG  .  7)  IR--3 
I F  ( IR . to . 0 ) I R  =  0 
IF  < IR . ED . 9) IR  =  1 
]  F  <  IR . EG .10)10-2 
IF  dR  - EQ .11) IR-3 
1 F <  IR  .EG.  12)100 
IFdR.EQ,  121IL-0 
IFdR.F.O, 12)10*0 
IFdC.NE . 1 >  GUTD  10 
It  dB.NE.0)  GOTO  10 
If (ANS.EG. 0 .0)  GOTO  25 

INITIALIZE  THE  LEX  I  DAI  A 

CALL  LIMIT 

CALL  DSMRG<0, 0,539. 639,0) 

CALL  DSCHNI-i , -1 , -1 > 

JUMP  TO  10  IF  INITIALIZATION  HAS  BEEN  ALREADY  COMPLETED 

CALL  DSL IM( 38+60* IR ,8+IL*7u ,69+60*IR ,39+IL*70  > 

CALL  DSPUTIIAl .1024) 

CALL  DSZOMU.1,2) 

CALL  DSSAO(31+60*IR ,45+IL*70 , lOOOB.O .1 > 

CALL  DSTXT  <  6 , IMD) 

CALL  DSSAO<  40+60*IR , 56+1 L *70 , 1 OOOP , 0 , i ) 

CALL  D3TXT  <3 , INUM) 

CONTINUE 

RETURN 

END 

THIS  SUBROUTINE  DISTINGUISHES  BETWEEN  A  BRIDCE  OVER  HATER  AND 
ANY  OTHER  TYPE  OF  GEOGRAPHY.  IF  A  BRIDGE  IS  DETECTED,  ANS-1  IF 
IT  IS  NOT.  ANS-0. 

SUBROUTINE  BRDGEl IA1 , ANS , ICOL , IROW) 

DIMENSION  TX<32 .32) , I A1 <  32 , 32 ) 

N-0 

ANS-0 

Nl-0 

ICOL-O 

IROH-O 

DO  10  1-1 ,32 
DO  10  J»1 ,32 
TX ( I  ,  J>«IA1 d  ,  J ) 

CONTINUE 
DO  20  J-1,32 
TX<1,J>»S 
TX(32 , J)*5 
TX(J,l>-5 
TX  <  J , 32 ) -5 


in 


o^vv 

,?0 

CONFIN'II 

1)600 

1X1  40  J=?  ,  31 

0601 

Dll  30  1=2.31 

060.? 

II  <  rX(l  ,J>HX(I-1  ,J>  ,Fl|.?>N=N+l 

060,5 

ir  <  TXIJ, I>*IX< J.I-1 ) .EU.?)N1=NI+1 

0604 

30 

CONI INOF 

060S 

IF  (N.Cr.  IH.AND.N.OF.  .?>H:OL«ICOL*l 

0606 

U (N1.L1 . IB. AND.Nl . GF . ? > I RUW^ IROU+ 1 

0607 

N^-0 

060R 

N1  =0 

0609 

40 

CtINI  INOF 

0610 

lr(ICm..RF.O)ANf5=l 

0611 

IF  ( IROU .  GF. .  8)  ANS=1 

061? 

RFIIIRN 

0613 

F.ND 

0614 

FND« 

